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The anisotropy and polarization of the cosmic microwave background radiation 
(CMBR) induced by the scalar and tensor metric perturbations are computed in the 
long-wavelength limit. It is found that the large-scale polarization of CMBR induced 
by the decaying tensor mode can reach a few percents. This is different from the 
scalar or inflation-induced tensor mode, whereas the polarization is at least two orders 
of magnitude lower. The effect of matter re- ionization on CMBR is also considered. 
We conclude that measuring the polarization of CMBR on large-angular scales can 
probe the ionization history of the early Universe, and test different cosmological 
models. 



PACS numbers: 98.70.Vc, 04.30.-hx, 98.80.Cq 



The recent detection of large-angular-scale temperature quadrupole anisotropy in the 
cosmic microwave background radiation (CMBR) by the DMR aboard the COBE satellite 
opens a window to our understanding of physics associated with the initial conditions of the 
early Universe. The most important conclusion to be drawn from the COBE data on the 
temperature correlation function is that they are consistent with a scale-invariant spectrum 
of primordial density perturbations predicted by inflationary cosmology 0. 

There are two sources of CMBR anisotropy, namely, the density perturbations (scalar 
mode) and the primordial gravitational waves (tensor mode). As density perturbations 
and gravity waves enter the horizon, they induce large-scale temperature anisotropics in the 
CMBR via the Sachs- Wolfe (SW) effect [^]. It is generally believed that the scalar mode can 
account for the CMBR anisotropy, which is regarded as evidence for the seeds of galaxies 
formation. 

Recently, it was argued that the anisotropy might be dominated by the tensor mode [^Q, 
and was showed that the tensor-mode dominance actually occurs in certain inflation models 
0. In Ref. 0, the authors suggested that by comparing large- and small-scale anisotropy 
measurements, one can separate the scalar- and tensor-mode contributions. At small-angular 
scales (~ 1°), the tensor-mode contribution to the CMBR anisotropy decreases relative to 
the scalar mode. A preliminary measurement of the CMBR anisotropy on 1° have indicated 
a low value when assuming a scalar- mode contribution p[. If the result should hold true, 
then this may be an evidence for the tensor-mode dominance. However, an alternative 
explanation is that late matter re-ionization may wash out the small-scale anisotropy too. 
A potential method for distinguishing tensor mode from re-ionization scenario is to measure 
the polarization component of the CMBR on small-angular scales: re-ionization leads to 
much greater polarization. Anisotropic radiation acquires linear polarization when it is 
scattered with free electrons (Thomson scattering) 0. The effect of Thomson scattering 
on the polarization of CMBR for the scalar-mode and tensor-mode cases can be found in 
Refs. [0 and [|lT| respectively. Ref. [10] showed that roughly 10% of CMBR anisotropy is 
polarized on angular scales less than 1°. Ref. [11] showed that tensor mode induces large 
polarization at small angular scales if primordial spectrum extends to small scales. In this 
paper we estimate the polarization to anisotropy ratio of CMBR due to long-wavelength 
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scalar and tensor modes in an universe with and without a late re-ionization era. The 
result shows that polarization measurements of CMBR on large-angular scales can probe 
the reionization history of the early Universe and possibly test different cosmological models. 

To study how polarized photons propagate in the Universe, one need to solve the equa- 
tion of transfer for photons. The general formalism for the subject of radiative transfer is 
given in Ref . . In general, arbitrarily polarized photons are characterized by four Stokes 
parameters, n = (n/, n,., n^, where n = rii + is the distribution function for photons 
with / and r denoting two directions at right angle to each other. 

We shall use the units c = h = 1 throughout. The metric that we use is of the fiat 
Robert son- Walker form 

ds^ = a\r]) [dr]^ - d^^) , (1) 

where 0(77) and drj = dt/a{t) are the scale factor and conformal time respectively. Then, it 
follows that, for a matter-dominated universe, 

p = 0. aM = ^. ,„ = 1. (2) 

where p, p, and Hq are the pressure, energy density, and present Hubble constant respectively. 
Note that we normalize the conformal time to unity today. In this metric, fitotai = ^dm + 
Qb = ^, where Qdm and Qb denote respectively the dark and baryonic matter. 

Then, the equation of transfer for an arbitrarily polarized photon is governed by the 
coUisional Boltzmann equation. 



^ d ^ d \ du dn 

^or] o:x. J drj ou 



n / / P{fi, (p, p , (j) )ndfi 



(3) 
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where cxr is the Thomson scattering cross section, Nf, is the number of free electrons per 
unit volume, (/i = cos6', 0) are the polar angles of the propagation direction e of the photon 
with a comoving frequency u, and P is the phase matrix for Thomson scattering given in 



Ref. \12 



The first term on the right-hand side of Eq. (3) describes how the frequency of a photon 
is shifted by metric fluctuations. The second term accounts for the Thomson scattering effect 
on photon polarizations. For an universe with the metric (1), the solutions of the perturbed 
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Einstein field equation are given in In the matter-dominated epoch, the solution for the 
spatial part of the metric perturbations is given by 

^ Wi (t) - - t) ^ ^ ^ ^ - To^" 

up to a gauge transformation, and the matter density fluctuation 

^"=32^5^ (if • 
where the functions A{x) and -B(x), Cj(x) and Dij{'x,ri) correspond to the scalar, vector 
and tensor metric perturbations respectively. In the presence of the metric perturbations, 
the rate of change of the photon frequency is given by the SW formula 

u di] 2 dr] drf 
At a sufficiently early epoch (onset of matter-dominated era at redshifts z < 10^, say), 
photon radiation is assumed to be a blackbody and unpolarized. Afterwards, the photons 
interact with matter through Thomson scatterings, and the subsequent evolution for photon 
distribution function is determined by Eq. (3). Solutions for Eq. (3) which are relevant to 
the present consideration have been obtained in two cases: the first case is that the frequency 
shift is owing to weak gravitational waves (corresponding to the Dij functions in Eq. (4)) 
TT| , and the second due to anisotropic expansion of an universe P,p!3[|. It was found that 



in both cases the polarization to anisotropy ratio of CMBR on large-angular scales is about 
a few percents. If matter re-ionization occurs subsequent to the hydrogen recombination, 
the ratio will be much enhanced in the anisotropic model ll3|. It is realized that the case 
of unequal Hubble expansion along three axes is equivalent to the infinite-wavelength limit 
of the tensor-mode case up to an overall normalization factor, which is determined by the 



magnitude of the gravitational wave |Tl|. Also, it is equivalent to an axisymmetric universe 
up to a spherical harmonic function |[T^. Below we will consider the effect of a generic 
scalar perturbation to the large-scale anisotropy and polarization of CMBR, i.e., setting 
Dij = Ci = (hence hjo = 0) in Eq. (4). Then, we will consider the large-scale anisotropy 
and polarization of CMBR due to the tensor perturbation (corresponding to Dij) generated 
in an inflationary cosmology. 

We assume a density perturbation wave 



5p = 5^p e^'^-^ (7) 

at the present time 770 = 1 with comoving wave vector k. It would induce a small amount of 
anisotropy and polarization in the CMBR. Thus we expand the photon distribution function 
as 

n = no + ^noni(7?) e*'^ , (8) 

where Uq is the blackbody distribution function, no = |no(l, 1,0,0) denotes unpolarized 
photons, and ni is a small fluctuation. Choosing k along with the z-axis, and making use 
of the axial symmetry of the problem, we obtain from Eq. (3) the equation of transfer for 
the reduced perturbation ni = {nn^nir) (ni^ = rii^ = 0) as 

\dr] ) dhiv V df] \\ 



— arN^a 



■ 3 .1 /2(l-/.2)(l_^'2^+^2^'2 ^2 

"^-8 III ^'^ J"^^^ 



(9) 



where 7 = (91nno/i91nz/ = —1 in the Rayleigh-Jcans region. 

To evaluate u'^du/dr] in Eq. (9), we separate two cases: A ^ and B — 0; B ^ and 
^4 = 0, corresponding to Sp/p which decreases and increases with time respectively. The 
increasing mode is believed to be responsible for the structure formation of the Universe. 
Prom Eqs. (4)-(7), we find that 

Idu f-3|f^/iV^-^ for^^O; 



^ . (10) 
udr] |^2^r7//V*-^ for S ^ , 

where po = SH^/ (SuG) is the present average energy density. In either case, Eq. (9) is 
too complicated to be solved. However, by taking infinite-wavelength limit, i.e., k 0, one 
could extract simple solutions from Eq. (9). These solutions are indeed good estimates of 
the anisotropy and polarization of CMBR on large-angular scales (see below). 

By substituting Eq. (10) in Eq. (9), taking A; — > limit (i.e., neglecting the term ikjini), 
and performing mode-decomposition, 

ni = «(/.^-M J+/^(1-/^')(_J ' (11) 



where a and (3 are respectively the anisotropy and polarization of the cosmic radiation at 
time 77 (in essence, a represents the temperature anisotropy, ST/T, and (3 represents the 
fractional temperature difference between two orthogonal directions along the hne of sight), 
we obtain 

| = -,AM-.,iV.a(^„+5^), (12) 
| = -.,iV.«(l« + H/3). (13) 

where, in terms of redshift z, 

2,2 1^(1 + ^)4 B^O. 

Integrating a and /3 in Eqs. (12) and (13) with respect to z, we obtain their present values 



a — —7 / A{z 



7 7 

1 1 



{l + z)~^dz (15) 



/3 = -7 / A(^) -e-^(") - -e-w-(-) (l + z)-idz (16) 
Jo 17 7 J 

where A{z) is given in Eq. (14) and the optical depth is given by 

^{^) =To [ Xe{z')\/l + z'dz' ; Xe = ^ , To = ^'^"^^ ^ (17) 
io JSIb nQUiB 

where Xe is the degree of ionization, rriB and Nb are respectively the proton mass and 
number density, pc is the critical energy density, and VLb = Pb/Pc- Henceforth, we assume 
Hq = 100/i km sec-^Mpc-i with h = I. Then, tq = 0.0696 VLbH. 

To evaluate a and (3 in Eqs. (15) and (16), we need to determine the optical depth which 
is determined by the dynamics of hydroden recombination and the history of re-ionization. 
Let us first assume an instantaneous recombination: Xe = 1 for z > Zr and Xe = for 2; < Zr, 
where z^ is the redshift at which hydrogen recombination occurs; z.,. = 1350, 1450, 1500 for 
flB = 0.01,0.1,1 respectively. Express a = ±q;o7 Sqp/pq, where the plus (minus) sign 
corresponds to the A {B 0) case. 

For A ^ we obtain from Eqs. (14)- (17), cto ~ 5x 10^ and /3/q; ~ 10"^- 10"^ for Qb from 
0.01 to 1. In this case, a ~ —0.5 when assuming a horizon-sized density perturbation with 
amplitude 5oP/Po — 10~^, which is too big when compared to observational data. However, 
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the A ^ case is equivalent to (i) an axisymmetric universe up to a numerical factor, and 
(ii) that of an infinitely long decaying gravitational wave up to a numerical factor times a 
spherical harmonic function. Nevertheless, all the three cases have the same (3/ a ratio. The 
reasons for the above equivalences are the following. For an axisymmetric universe with the 
present shear given by AHo, the rate of change of the photon frequency is given by []T^ 



Idu _ 2AHo 1 1 ^ 

u dr] Hq 7]^ 3 ' 



(18) 



which is similar to the A case in Eq. (10) {k —>■ 0) up to an irrelevant monopole term. 
In the presence of a monochromatic gravitational wave with comoving wave vector k and 



polarization A = + mode, it was shown that fTT 



Idu _1 
udr] ~ 2^ 



/i ) cos^ 



-ikx 



kr]. 



(19) 



drj ' ' 7]^ 

where / is a constant wave amplitude. Taking the — limit, Eq. (19) is similar to the 
A ^ case in Eq. (10). For such wave, the mode decomposition is given by 



ni 



-(1 
2^ 



/i^)cos20 
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1 













(20) 



( (l + /i2)cos20 \ 
;i+;u2)cos20 . 
4/isin20 / 

For 5 ^ 0, we find ao - 1 and (3/a ^ 10"^ - 10"^ for VLb from 0.01 to 1. Then, 
a ~ Sqp/po — 10^^ which is the well-known SW contribution to the large-scale anisotropy 
of CMBR. jS is below 10"^" which is consistent with the values of lower multipole moments 



that were given in Ref. |10 



To determine the effects of non-instantaneous recombination, we approximate the history 
of ionization by, 

2 < 900 , 

(21) 

1 z > 1500 , 



Xe 



3.2 X IQ-^rij^^ 
5.68 X 10^Q]^^z~^e~^ 900 < 2 < 1500 



where we assume that a residual ionization remains for z < 900. For 900 < z < 1500, we 



use the results given in Ref. [14] which take into account the possibility of two-photon decay 
of the hydrogen 2S level during the process of recombination. For z > 1500, we assume a 
complete ionization (an exact treatment would require using the Saha equation). Plugging 



Eq. (21) in Eq. (17) and then using Eqs. (14)-(16), we find that the degree of anisotropy, on 
comparing with the instantaneous case, remains approximately the same, whereas the degree 
of polarization can possibly greatly enhanced. Being weakly dependent on the value of Qb, 
the ratio (3 /a increases to about —4% for the case of A ^ and — 10~^ for that of -B 7^ 0. 
These results are shown in Figs. 1 and 2. Our value of the degree of polarization for the 
A case is consistent with the 2% result that was given in Ref. [0. To understand our 
results for the non-instantaneous case, we examine the value of the optical depth evaluated 
at the beginning of hydrogen recombination. This value denotes the number of Thomson 
scatterings between epochs z = z^. and z = 0. We find that t{z = 1500) ~ 90, which is 
much larger than one. This explains why the polarization is sensitive to the recombination 
dynamics. 

Finally, we consider the effects of matter re-ionization subsequent to the recombination 
epoch. This re-ionization of matter may attribute to the reheating by radiation released from 
protogalaxies during galaxy formation. To demonstrate this, we approximate the ionization 
fraction hj Xe = ^ for z > Zr] for Zh < z < Zr] I for z < Zh, where Zh is the redshift at which 
the re-ionization occurs. Note that the case of instantaneous recombination is equivalent 
to setting z^ = 0. Using Eqs. (14)-(17) and varying Zh from 10 to 100, we find that the 
polarization, being very sensitive to the re-ionization, becomes much larger (see Fig. 2). To 
understand why the polarization curves all suddenly change around Zh = 10, we evaluate 
the optical depth ai z = 10. We find that t(10) ~ 1.6^2^. This explains the jumps for all 
the cases except the B ^ and ^Ib = 0.01 case. However, in this case the polarization 
before recombination is extremely small, even a small amount of scatterings during the re- 
ionization epoch can induce an increase in the polarization over its prerecombination value. 
As for the anisotopy Oq, it is found that it remains almost unchanged in the case of B ^ 
and it is diminished in that of A ^ (see Fig. 1). As mentioned above, the A ^ case 
is equivalent to an axisymmetric universe, this decreasing behavior agrees with the results 
given in Ref. ||T3| . 

We are now to justify the infinite- wavelength approximation in our calculations. Since 
horizon-sized waves dominate the contributions to the anisotropy and hence polarization of 
CMBR on large-angular scales via the SW effect either for the scalar ||T^ or tensor mode 
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ip. Thus, one should compare the Thomson scattering mean-free path, {aTNea)~^ , with the 
wavelength of horizon-sized density perturbations, kf^l^.^ at the moment when the contribu- 
tions to the integrals in Eqs. (15) and (16) are dominant. If the former is much smaller than 
the latter, then the term ikfini in Eq. (9) can be neglected. For example, in the case of in- 
stantaneous recombination, the polarization j3 receives contribution mainly from z > Zr- For 
Qb = 0.01 and Zr = 1350, axNea > 2500 for z > z^ which is much bigger than k^or — 1- For 
other values of Qb and Zh, we find that axN^a is at least comparable to or much bigger than 
1. Thus, we conclude that taking the infinite- wavelength limit is an efficient approximation 
in calculating the lower multipole moments, especially the quadrupole moment. However, 
the finite-wavelength correction is of interest for intermediate-scale calculation. 

We now turn to the infiation-induced tensor mode calculation. The calculation is similar 
to the case of decaying tensor mode. The rate of change of the photon frequency is given 
by Eq.(19) with he^''^ replaced by 3A{k)ji{kr])/{kr]), where A^lk) = 8i;/(37r) is the scale 
invariant spectrum, ji{kr]) is the spherical Bessel function of order 1, and f is a parameter 
depending on the infiation scale By taking the infinite-wavelength limit, the rate of 

change of the photon frequency decreases to zero. This is expected since the amphtude of 
the tensor mode well outside the horizon is constant. Therefore, the above limiting approach 
fails in this case. However, as we have mentioned that horizon-sized waves dominate the 
contributions, we instead proceed the calculation by simply dropping the term ikfiiii and 
taking k = 2tt in Eq. (9), and using t> = 4 x 10^^^. The results are shown in Figs. 1 
and 2. The behavior of the results is very similar to that of the B ^ case. This makes it 
difficult to distinguish the scalar mode from the infiation-induced tensor mode by measuring 
the large-scale anisotropy and polarization of CMBR. We have also carried a full numerical 
calculation of Eq. (9) by taking ~ 1 and found that dropping the term ikfiiii does not 
make any significant difference. 

In conclusion, we have calculated the large-scale anisotropy and polarization of CMBR 
due to the scalar and tensor metric perturbations, as well as the effects of the non- 
instantaneous recombination and matter re- ionization. The results are shown in Figs.l 
and 2. The ao (up to an overall normalization factor which is determined by the amplitude 
of the decaying tensor mode) and (3 /a values for the A ^ case also represent respectively 
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the anisotropy and the polarization to anisotropy ratio due to the decaying tensor mode. Al- 
though it seems unlikely that measurable super-horizon sized decaying tensor modes could 
be generated in the early Universe, however, the large polarization due to this decaying- 
type perturbations or anisotropic Hubble expansion |16| remains an interesting possibility. 
It is shown that measuring the degree of polarization of CMBR on large-angular scales is a 
potential method for probing the ionization history of the early Universe. The linear polar- 
ization of CMBR was measured at 33 GHz over 11 declinations from 37''S to 63°N. Fitting 
the data to an axisymmetric model and spherical harmonics through the third order yields 
\(3\ < 3 X 10~^ and |/?| < 6 x 10~^ respectively at 95% confidence level Jl^. These upper 
limits can hardly be used to put any constraint on the present calculations. Currently, the 
linear antenna aboard COBE has collected data for the polarization component of CMBR 
which is being analysed. Any trace of polarization or a better limit will prove invaluable to 
our understanding of the early Universe. 
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FIGURE CAPTIONS 



Figure 1. Large-scale anisotropy of CMBR versus the redshift at which matter re- ionization 
occurs. Symbols A, B and T denote respectively the cases of A 7^ 0, S 7^ and the 
inflation-induced tensor. For the cases of A 7^ and B 0, a — ±q;o7 5op/po. For the 
tensor case, a = 700 10^^. Zh = corresponds to instantaneous recombination. Each solid 
triangle corresponds to non-instantaneous recombination for Qb = 1,0.1 and 0.01. Dashed, 
solid and dotted curves correspond respectively to fls — 1, 0.1 and 0.01. 

Figure 2. Large-scale polarization to anisotropy ratio of CMBR. Symbols A, B and T 
denote respectively the cases of A 7^ 0, S 7^ and the inflation-induced tensor. Zh — 
corresponds to instantaneous recombination. Each solid triangle corresponds to non- 
instantaneous recombination for CIb = 1,0.1 and 0.01. Dashed, sohd and dotted curves 
correspond respectively to Hb = 1,0.1 and 0.01. 
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